Iterations of nonlinear entanglement witnesses 
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We describe a generic way to improve a given linear entanglement witness by a quadratic, nonlinear 
term. This method can be iterated, leading to a whole sequence of nonlinear witnesses, which become 
stronger in each step of the iteration. We show how to optimize this iteration with respect to a 
given state, and prove that in the limit of the iteration the nonlinear witness detects all states that 
can be detected by the positive map corresponding to the original linear witness. 
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I. INTRODUCTION 



When Erwin Schrodingcr introduced the notion of en- 
tanglement for certain bipartite quantum states in the 
thirties of the last century, he might not have imagined 
that nowadays this peculiar phenomenon constitutes the 
fundamental resource for such fascinating tasks like quan- 
tum cryptography or quantum teleportation. By defini- 
tion, an entangled state is not separable, which means 
that it cannot be prepared by local operations and clas- 
sical communication Any separable state p^]^ can 
be written as the convex combination of pure product 
states, i.e.. 



(1) 



with a probability distribution {pi\i and corresponding 
pure states l^/'i)^ and for the local subsystems [41 



Despite its importance for the field of quantum infor- 
mation theory, the properties of entangled states are not 
fully explored yet. Even to determine whether a given 
quantum states is entangled or not is still an open prob- 
lem, although considerable progress has been achieved 
along this directions over the last decade, see Ref. 0. In 
fact this so-called separability problem can already be re- 
garded as a research field on its own and several different 
results have provided insight into this problem: the for- 
mulation of operational criteria which are sufficient to de- 
tect either entangled or separable states d, H, H, d, 0j 111 j 
different ways to tackle the problem by numerical means 
[13, [HI , or the reformulation of the separability prob- 
lem into a different contex t [13 . [isj . Remarkably, in the 
case of low dimensionality , or for a p articular class of 
even infinite dimensional states [l3j[lM, separability 
problem is solved. 

Another approach to the separability problem inves- 
tigates entanglement witnesses [1, [l^ [l3l- An entan- 
glement witness W is an observable which has a non- 
negative expectation value on all separable states, and 




FIG. 1: Because the set of separable states (dotted) forms 
a convex subset to the set of all possible density operators 
(outer, thick-lined), any entangled state can be separated by 
a corresponding linear entanglement witness Wi . The present 
manuscript describes an iteration method to successively im- 
prove entanglement witness, corresponding to a certain posi- 
tive but not completely positive map, by quadratic nonlinear 
terms, and thus generates a sequence of nonlinear entangle- 
ment witnesses W2, W3, . . . VK> of increasing strength. 



therefore any negative expectation value signals the pres- 
ence of entanglement. Those kind of operators offer a 
powerful tool to verify the creation of an entangled state 
in an actual experiment, since one only has to measure 
the corresponding observable, cf. Refs. [isl. [lOl. [20j . 

Witnesses allow the following geometrical interpreta- 
tion: As schematically shown in Fig. [T] the set of all pos- 
sible bipartite quantum states forms a convex set, while 
the restricted set of only separable states constitutes a 
convex subset of it. Therefore it is possible to separate 
any given entangled state from the set of separable states 
by a corresponding hyperplane. This hyperplane repre- 
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sents the corresponding entanglement witness, or, more 
precisely, the states for which the mean value of the wit- 
ness equals zero. The witness detects the entangled state, 
because its expectation value allows distinction between 
the states from either the "left" or the "right" hand side 
of the witness. 

From this geometrical picture, it is natural to ask 
whether it is possible to improve a given linear entan- 
glement witness by some nonlinear correction in order to 
approximate the set of separable states better. Indeed, 
as shown in Refs. [2ll,[22| this is always possible. In these 
references a general method has been provided to com- 
pute a nonlinear improvement to a given linear witness. 
The nonlinear corrections consist of terms which are typi- 
cally squares of certain expectation values, and the whole 
expression is still positive on all separable states. The 
nonlinear witnesses may also be viewed as linear observ- 
ables, but acting on several copies of the input state [23l |. 

In this paper we provide an extension of this improve- 
ment idea which leads to a whole sequence of nonlinear 
witnesses. We first present a method to compute a non- 
linear improvement for a given witness. This method 
can then be iterated, leading to a sequence of nonlin- 
ear witnesses. In this iteration, the nonlinear witnesses 
contain higher nonlinearities and become stronger in each 
step. Each iteration step requires a choice of an arbitrary 
unitary transformation. We show how to optimize this 
choice with respect to a given state, in the sense that this 
state should be detected after as few iterations as possible 
and show that this iteration finally detects all the states 
which are detected by the positive map corresponding to 
the witness. 

In detail, our paper is organized as follows. In Sec. |TT] 
we present the main idea of the iteration and fix our no- 
tation. We consider witnesses from the criterion of pos- 
itivity of the partial transpose (PPT) Q at this point, 
however, as we will see at the end, the methods can di- 
rectly be used for arbitrary witnesses using the duality 
between witnesses and maps [13, [H, |2^. In Sec. IIIII we 
discuss the optimization of the iteration with respect to a 
given state and prove its main properties. In Sec. IIVI we 
discuss what happens if the iteration is not optimized and 
the unitary transformation is chosen randomly. Interest- 
ingly, it turns out that the resulting nonlinear witnesses 
are still strong. In Sec. |V] we present several examples 
of our iteration. Finally, we discuss extensions of our 
method to arbitrary witnesses, and conclude. 



II. ITERATION PROCESS 

The main idea of the iteration process is a particu- 
lar way of improving a given entanglement witness by 
some quadratic, nonlinear term. The point is that this 
improvement method has the advantage that it can be 
employed several times. This leads to a sequence of en- 
tanglement witnesses of increasing strength. 

In this section we will first introduce the basic idea of 



the iteration, and then explain some generic properties 
of the iteration process. By giving a simple example se- 
quence, we can show that the iteration process allows to 
"curve" a given entanglement witness in such a way that 
it will detect an arbitrary given state violating the PPT 
condition at some point in this sequence. 



A. Main idea 

Let us explain the main idea in the simple setting of 
witnesses coming from the PPT criterion. For a given 
density matrix p = u PijM |u) {kl\ the partial trans- 
pose with respect to the second system is given by 

^J2P^^^k,\^J){kl\. (2) 

The celebrated PPT criterion states that if p is separable, 
then its partial transpose is positive scmidefinite, i.e., it 
has no negative eigenvalues. Hence any quantum state 
with a negative partial transpose (NPT) must necessarily 
be entangled. 

This criterion has two consequences. First, for a sepa- 
rable state we have ((/)| p^ |(/)) > for any vector \(j>) . This 
is equivalent to ti{p^ P) > for any positive semidefinite 
operator P or to ti{p^ BB^^) > for any matrix B. (See 
also [131. ) Clearly, these conditions are nothing but a 
reformulation of positive semidcfinitcness of the partially 
transposed state. 

Second, the PPT criterion allows to construct directly 
a witness for a given state violating it. Namely, if p^ 
has a negative eigenvalue A_ < and \tp) denotes the 
corresponding eigenvector, then 

Wo = m {^Pf (3) 

is a witness detecting p, as tr(pVFo) = tr{p\ijj) (tpf ) ~ 
tr{p^ \ip) = A_, while for all PPT states, thus also 
separable states, it holds tr(pVFo) > 0. Note that more 
generally any observable of the type W = [BB^)^ is posi- 
tive on all PPT quantum states, and any NPT entangled 
state can be detected by a witness of this form. Our 
aim is to improve this witness by a sequence of nonlinear 
witnesses. 

To start, note that B might be a linear combination 
of some other operators {-Fiji, i.e., B = J^i'^i^i- This 
implies that for separable p one has 

trip^BB^) = ^'(P^P'^Pj) ^ 0- (4) 

As B is not fixed yet, this has to hold for all possi- 
ble combinations of ci € (D. Such a condition can be 
fulfilled if and only if the matrix M{p^) defined above 
is positive semidefinite itself. This particular matrix, 
which we call expectation value matrix (EVM) in the 
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following, has recently drawn much attention in differ- 
ent contexts of the quantum information literature, see 
Refs. [ll,[23,[2i,[2i,[33,[3l|. 

In order to improve the linear entanglement witness, 
we consider the EVM for the specific operator set Fi = 
and F2 = B, where U constitutes an arbitrary uni- 
tary and the operator B characterizes the linear witness 
W = {BB'^)^. Assuming normalization of the partially 
transposed state tr{p^) = 1, the corresponding EVM be- 
comes 

M(o^)-( ^ tr(pr(sc/)t) \ 
^^^^P ^ ~ [trip^BU) trip^BB^) J' 

As long as p is separable, p^ > forms a positive 
semidefinite operator, and hence the corresponding EVM 
M{p^) > will be positive semidefinite as well. To fi- 
nally achieve the quadratic, nonlinear improvement, one 
employs the Sylvester criterion (s^ . which states that 
a given matrix can only be positive semidefinite if all 
principle minors are non-negative. This results in condi- 
tions on the diagonal elements and the determinant of the 
EVM given by Eq. O Hence the EVM for any separable 
state must fulfill 

Wi^{p) = trip^BB^) > 0, (6) 
M^nl(p) = tr(/BSt) - \tT{p^BU)\' > 0, (7) 

since otherwise the partially transposed state has a least 
one negative eigenvalue and therefore signals the presence 
of entanglement of the corresponding quantum state. 

El- 

The first condition given by Eq. [6] represents the origi- 
nal linear witness Wl- The nonlinear entanglement wit- 
ness Wnl defined by Eq. [7] is a strict improvement of 
the linear witness, as one further subtracts another non- 
negative term [43 |. 

The main point of this construction is that the nonlin- 
ear improvement Wnl can again be written as an expec- 
tation value of a linear witness, and hence the procedure 
can be iterated. Namely, if we define 

B' = B'{B, U, /) := BU - tr(p^BC/)l, (8) 

we have 

W^nl(p) =tr(pi^B'B't) (9) 

and we can write down a further improvement and so 
forth. Note that the defined operator B' depends on 
the starting operator B, the chosen unitary U and the 
considered partially transposed state p^ . 

To give an example, for the case of the witness Wq 
in Eq. [3] we can choose B = [ip) {ipl, and an arbitrary 
10) = [/t If/)), then the nonlinear entanglement witness 
becomes 

W^nl(p) = tr{pWo) - I tr [pi\^P) (cP\f] \\ (10) 

which is similar to the quadratic improvement from 
Ref. [H. (See also [H.) 



This iteration poses now a set of interesting questions. 
First, one may ask whether one can always find a se- 
quence of improvements of a witness W, which detects a 
given state p, which is not yet detected by the starting 
witness W. As we will see, this is the case. Second, one 
may ask how one can find the optimal iteration, such that 
the nonlinear term which is subtracted is maximized in 
each step. Third, one can investigate what happens if 
the unitaries are chosen randomly. We will address all 
these questions in the following. 

Finally it should be noted that the described iteration 
process to obtain nonlinear witnesses only depends on the 
knowledge whether the partial transposed state is posi- 
tive semidefinite or not. Therefore no entangled quantum 
state with a positive partial transpose is detected by any 
of those nonlinear entanglement witnesses, however the 
method can directly be generalized to other positive but 
not completely positive maps, cf. Sec. IVII 

B. Definition, properties and example iteration 

The following definition formally introduces the itera- 
tion process. More specifically the iteration contains the 
sequence of entanglement witnesses, which arc defined 
via their relation to the operators B (cf. Sec. Ill Ap . the 
corresponding sequence of expectation values and the re- 
sulting sequence of improvements. However, note that 
the iteration process itself is not fully specified yet, be- 
cause so far no explicit recipe to draw the unitary oper- 
ators has been fixed. In general, one could distinguish 
different kinds of iterations depending on the strategy 
to choose each unitary. In the subsequent sections we 
mainly distinguish three different iterations: the opti- 
mized iteration (Sec. IIIip in which each unitary is cho- 
sen such that it maximizes the improvement for a given 
quantum state, a random iteration where each unitary is 
chosen at random and the averaged iteration (Sec. IIV|) . 
which considers the averaged expectation value taken 
over all possible unitaries. 

Definition II. 1 (Iteration process). Given the operator 
p^ with ti{p^) ~ 1 together with a starting operator Bi, 
the iteration process contains the following sequences with 
n > 1. Each sequence depends on the set of unitaries 
Ui, . . . Un-i, which are chosen in the steps before. 

• The sequence of operators {-B„}„, which are recur- 
sively defined by 

Bn = Bn{Bn-l,Un-l, P^) 

B„_i[/„_i-tr(pi^B„_i;7„_i)l, (11) 
and the operator Bi is the given starting operator. 

• The sequence of expectation values {wn}n , that are 
defined via 

Wn = Wnip"^, Bn) := tr(pi^S„Bt ). (12) 
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• The sequence of improvements {c„}„, which char- 
acterize the detection improvement in each step, are 
given by 



Cnip^ ,Bn,Un) := \tT{p^ BnUn)\ 



(13) 



Note that each operator Bn = Bn{p^ ,Ui, . . .Un-i), 
each expectation value Wn = Wn{p,Ui, . . .Un-i) and 
each improvement c„ = Cn{p^ , Ui, . . . , Un-i) depends on 
the unitaries J7i,...J7„-i chosen in the iteration steps 
before, the operator p^ and the given starting operator 
Bi] however, this exphcit dependence is suppressed in 
most of the following cases. Even without specifying the 
exact sequence of unitary operators, the iteration process 
has some generic properties: 

Proposition II. 1 (Generic properties). For the iteration 
process, introduced by Def. \II.1[ one has the following 
generic properties, which are independent of the starting 
operator Bi : 

• For a given operator p^ that satisfies tr{p^) = 1, 
the sequence of expectation values {wn}n is mono- 
tonically decreasing. 

• // the state is PPT, p^ > 0, then the expectation 
values Wn{p^) > are non-negative, for all possible 
n > 1. 

Proof. Using the recursion formula for Bn and the defini- 
tion of the c„, one gets an analogous recursion statement 



(14) 



The first statement directly follows from that. The sec- 
ond point follows from the definition of the expectation 
value Wn = tT{p^ BnB}^), and the fact that BnB}^ > 
and > 0. □ 

Before considering more specialized situations, we dis- 
cuss a simple example iteration. It shows that if one 
starts with an arbitrary entanglement witness corre- 
sponding to the operator Bi > it is always possible 
to improve the given witness in such a way that it de- 
tects the preselected NPT entangled state. 

Proposition II. 2 (Example sequence). For a given de- 
tectable operator p'" ^ with tr{p^) = 1 and Bi > 0, 
consider the following sequence of unitaries. 



Ui=P+-P-, Un 



-P+ + P_ V??. > 2, 



(15) 



where P_ denotes the projector onto the subspace of neg- 
ative eigenvalues, and P+ its orthogonal complement. 
Then there exists an Nq g N_|_ such that Wn < for all 
n > No, hence the corresponding sequence of expectation 
values {wn}n will detect the state; in fact the sequence 
{wn}n diverges to —oo. 



Proof. The recursion formula for the sequence of expec- 
tation values {wn}n, given by Eq. [M] allows us to prove 
the proposition by showing that the sequence of improve- 
ments {cn}n does not converge to zero, i.e., lim„^oo c„ > 
0; in fact, we even prove divergence of this sequence. 
Because of the special set of unitaries, each improve- 
ment is given by c„ = [tr(p'" £?„[/„)] . The first two 
improvements can be directly computed, and are given 
by ci = [tr{\p^\Bi)]'^ > and 



C2 = [tv{\p^\B,)\\p^\\,-tr{p^B,f 
> \\p^\\l{tr[{\p^\-p^)B^]}'>Q, 



(16) 



where ||...||i denotes the trace norm. The first inequality 
is valid for the case that tr(p'"i?i) > 0; if this is not the 
case, then C2 > trivially holds. The strict inequality 
in Eq. [16] comes from the fact that {\p^\ — p^) > and 
the assumption that Bi > 0. Note that in the case of a 
positive semidefinite operator p^ = \p^\ > 0, hence the 
second improvement vanishes. In the remaining part of 
the proof, we want to show by induction that 



II r||2(n-2) 
C2\\P 111 % 



V??. > 2. 



(17) 



The starting case n = 2 is trivial, and we only need to 
care about the induction step ni-^ n + 1. One arrives at 

C„+i = [tr{p^Bn)-tT{p^Un+l)tr{p^BnUn)f 

= {tr(pi^P„_iC/„_i)[l-tr(pr)] 
+ 11^^111 tr(p^P„C/„)}' 



II r||2[(n+l)-2] 
C2IIP 111 



(18) 



where we used the recursion formula for the operators 
Bn+i and UnUn+i = 1 in the first line. In the second 
step, we employed the recursion formula for the operator 



Bn and the identity tx^p^Un+i) 



while in the 



last step we used tr(p ) = 1 and the induction hypothesis 
in the end. Since Hp'^Hi > 1, the corresponding sequence 
of improvements {c„}„ diverges and the proposition is 
proved. □ 

This sequence of unitaries can always be "started" at 
each step of an arbitrary sequence; detection and diver- 
gence are ensured as long as the second expectation value, 
after one has used this particular sequence, is strictly 
greater than zero. 



III. OPTIMIZED ITERATION 

Suppose that one starts with a certain entanglement 
witness and that one likes to verify the entanglement 
of a given target state. Then it is of course desirable 
to detect the state with the least number of iterations. 
One possible way to achieve this goal might be to try 
to optimize the improvement for the specific quantum 
state in each step. Such a particular iteration method 
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is termed optimized iteration and is formally introduced 
in the following definition. Is is important to note that 
the sequence of unitaries, which are chosen during the 
iteration process, is mainly determined by the given tar- 
get state. From a geometrical point of view, this process 
corresponds to the task that one wants to curve a given 
entanglement witness along a certain direction, and the 
direction is linked to the quantum state. The final re- 
sult of this section is the statement that this optimized 
sequence will finally detect the initial target state. 

Definition III.l (Optimized iteration). For a given op- 
erator with tr(p'") = 1 and a certain starting operator 
Bi, the optimized iteration is defined as the iteration in 
which each improvement is maximized over the chosen 
unitary. The maximization of Cn — | tr(p'"i3„J7„)p can 
be obtained as follows: For any valid singular value de- 
composition p^Bn = VnDnW^ with £)„ > 0, one selects 
the optimal unitary 11°^^ = WnV^l, resulting in 

c„ = [tr{p^Br.U:^'f ^ [tr(AO]' = Wp^BnWj. (19) 

This particular unitary has the additional feature that 
p^ BnU°^*^ > is a positive semidefinite operator fjdi l. 

The extra requirement that the operator after the op- 
timization p^ BnU°^^ > forms a positive semidefinite op- 
erator is of particular importance in order to show detec- 
tion of an arbitrary target state, because it allows some 
predictions of the structure of the next optimal unitary. 
The proof itself relies on the idea that the sequence of 
improvements does not converge to zero. 

Theorem III.l (Detection for the optimized iteration). 

For any detectable state p^ ^ with tr(p'") = 1, and any 
strictly positive starting operator Bi > 0, the optimized 
iteration process will always detect the state. 

Proof. Because of the recursion formula of the sequence 
of expectation values, given by Eq. 1141 it suffices 
to prove that each improvement c„ is bounded from 
below by a strictly positive constant. For the op- 
timal iteration method, each improvement simplifies 
to c„ = dl with dn = tr(p^B„J7°P*) e M, hence 
it suffices to prove this bound only for the sequence 
{dn}n- Note that the first improvement, given by 
dl = tr(pi'Bif7°P*) = ll/Billi > 0, is strictly positive. 
According to the definition of the optimized iteration pro- 
cess, the optimal unitary in the first step U°^^ is deter- 
mined by a singular value decomposition of the operator 
p^ Bi and that one obtains a positive semidefinite oper- 
ator p^BiU°^* > 0. In the following, we want to prove 
the bound 

d„ = tr(/S„[/°P') > ciA„,i„ > 0, Vn>2, (20) 

where Amin denotes the absolute value of the most neg- 
ative eigenvalue of the operator p^ . The corresponding 
eigenvector is labeled by 

The case n = 2 already covers the main idea. Be- 
cause of the particular unitary that we chose in the first 



step, the next operator p^ B2 = p^BiU°^^ — dip^ is al- 
ready hermitian. Therefore we can expand the opera- 
tor in terms of its spectral decomposition as p^ B2 = 
12 Vi l^j) {"^il- From this explicit decomposition one can 
directly infer the next optimal unitary, which is given by 
^2^* — X]^iS^^('ft) ('^d; ^^'^ sign(77i) denotes the sign 
of the corresponding eigenvalue. As a result, one ob- 
tains a new positive semidefinite operator p'"_B2t^2^* = 
S \Vi\ > 0. These properties allow the following 

statement, 

d2 - trip^B2U°^') > {tlj\p^B2U°^'\ij) 

= J2 l^^l ^ E'^' 

> rflAmin. (21) 

In the first step, one uses that the trace over a positive 
semidefinite operator is lower bounded by the overlap 
over only one possible state. Note that the inequality 
(^ijj\p^ B 2112''^^ lip) > (^ijj\p^ B2\ip) only relies on the partic- 
ular spectral decomposition of the operator p^ B2 and the 
explicit choice of the unitary J/j ^* . Therefore this prop- 
erty holds in any step of the optimized iteration process. 
Using this idea multiple times enables us to prove the 
general statement, 

d„ = tr(/B„C/°P*) > (V^|p^B„-ii7°P\|V^) + d„_iA„,i„ 

> {^\p^ B„^2U°%\^^) + {dn~2 + d„-l)A„,in 

> . . . > (^|pi^BiC/°P'|7/-) + A,„i„(rfi + . . . + d„_i) 

> dlA,nin. (22) 

This finally proves the bound of Eq. [20l and therefore 
shows that the general sequence of improvements {c„} 
does not converge to zero. As a consequence the cor- 
responding sequence of expectation values {wn} will di- 
verge to —00, so that there exists a particular A^o G N_|_, 
such that Wn < for all n > Nq. Therefore the state 
^ will be detected at some point in the sequence. □ 

This theorem can be extended to the case in which 
the operator Bi > is positive semidefinite and ci = 
llp'^iJilli > 0, since the method to lower bound each im- 
provement does not rely on strict positivity of the start- 
ing operator. However in the specific case of ci =0, the 
partially transposed state p^ and the starting operator 
Bi are acting on completely orthogonal subspaces, and 
hence the sequence can never detect the state. 

IV. AVERAGED ITERATION 

Obviously, one drawback of the optimized iteration 
process is that the target state has to be known in ad- 
vance; however in certain cases such prior knowledge 
might be unavailable. Therefore, it is interesting what 
happens if one chooses the unitaries in the iteration pro- 
cess in a different, state independent way. A first simple 
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method would be to choose each unitary at random. At 
first sight, such a random iteration seems to produce only 
small improvements, however if one repeats the iteration 
many times one still can hope to detect many states. We 
will discuss this random iteration process with an exam- 
ples in the next section. 

In the present section we study the averaged iteration, 
in which, instead of using only a single sequence of uni- 
taries, one takes the average over all possible unitaries in 
the iteration process. As before, a negative expectation 
value for these mean values signals the presence of entan- 
glement. However, one should mention that the resulting 
nonlinear entanglement witnesses can not be written as 
a single nonlinear witness in the original form of the it- 
eration. 

The final theorem of this section states that there is a 
sequence of nonlinear entanglement witnesses, which in 
the limit of infinite many improvements, detects all NPT 
entangled states at once. Although this results seems 
surprising at first, such a method is already known 
and relics on the seminal spectrum estimation method 
introduced in Ref. |34| . Nevertheless, a similar statement 



can be derived using the idea of the averaged iteration 
process. 

Definition IV. 1 (Averaged iteration). For a given oper- 
ator that acts on an d- dimensional, composite Hilbert 
space, with tr(p'") = 1, and for a certain starting operator 
Bi, the averaged iteration defines the special sequence of 
expectation values {mJ„}„, in which one takes the aver- 
age over all possible unitaries Ui, . . . , Un-i, that one can 
choose up to the n-th step. More precisely, one defines 



Wr. 



dUi... dUn-iw„{p^ , Bi, C/i, . . . , C/„_i), (23) 



with Wn{p^ , Bi,Ui, . . . , Un-i) being the expectation value 
defined by Eq. WA and where dUi, for i = 1, . . .n, repre- 
sents the Haar measure of the unitary group U (d) with 
normalization J dUi = 1 . 

The intuition behind the detection statement of the 
averaged iteration process is very simple: One knows 
already certain sequences of unitaries {e.g. example or 
optimized iteration) that will detect a given quantum 
state and even diverge to — cx). In contrast for any other 
sequence of unitaries, in particular those which never 
detect the state, the corresponding expectation values 
are monotonically decreasing and bounded from above. 
Hence if one performs the average over all possible se- 
quences of unitaries, the resulting sequence of mean val- 
ues should diverge as well. 

Detection for the averaged iteration process is shown 
in a two step procedure: In a first step one obtains a re- 
cursion formula for the sequence of averaged expectation 
values. This particular formula simplifies the final proof 
of detection. 

Proposition IV. 1 (Recursion formula for the averaged 
iteration). In the averaged iteration process, each expec- 
tation value can be written as Wn{p^ , Bi) = tv{AnBiB\), 



for the starting operator Bi. The sequence of operators 
{An}n: that Only depend on the given input operator p^ 
with tr{p^) ~ 1, is recursively defined by 



Arn 



A„(/ 



An-1 



+ - [tr(A„_i)(/)2-(pr^„_i + A„_ipi^)],(24) 



starting with Ai = p^ , and d denotes the dimension of 
the underlying Hilbert space. 

Proof. The proposition follows by direct calculation of 
the averaged expectation values. In order to perform 
this task, we employ the identity that 



J dU dtr{AU) iv{BU^) tv{AB), 



(25) 



holds for arbitrary operators A, B. This can be proven for 
example by the Peter- Weyl theorem, see e.g. Ref. [ssj . 
(See also note [i3l-) To actually compute the averaged 
expectation value Wn, given by Eq. [23l one uses the 
following strategy: 

One starts with the average over the last chosen uni- 
tary U„-i. Using the recursion formula for the operator 
B„, and the given identity, it is possible to express the 
next integrand as the expectation value of a witness in 
the step before with a new operator A2 = A2{p^), hence 
/ dUn~itr{p^ BnB}^) = tr(A2-B„_ii?^_;^). One continues 
with the integration over the next unitary Un-2, and ex- 
ploits the same trick again to obtain a new operator ^3. 
In this step, one needs to be careful since the operator 
A2 might not be normalized any more tr(yl2) 1. In the 
end, one uses this idea exactly (n — 1) times until one 
has performed all the integration and ends up with the 
final operator An. For this operator one needs to com- 
pute the expectation value with the first witness BiBI, 

and obtains the final result Wn = tr(A„i3ii?j). 

Now that the strategy is fixed, one is left to perform the 
integration over only one unitary in order to obtain the 
recursion formula. Suppose the integrand is tr(Ai3„i?j^), 
where A denotes an arbitrary hermitian operator. Us- 
ing the recursion formula for the operator given by 
Eq. 1111 the integrand is expanded into 

tr(AS„Bi) = tiiABn-iBl_,) 

+ tr(A)|tr(/B„_iC/„_i)|2 

- tiiABn-iUn-i) tr{p^Ul_,Bi_,) 
-c.c. (26) 

Since the operator Bn-i is independent of the last uni- 
tary Un-i, one can directly perform the average over this 
last unitary. Each term can be integrated separately, and 
by using the given identity, this results in 

J dUn-i triABnBl) (27) 
= tr(^|A + i [triA)ip^r-{p^A + V)]|i?n-iBLi) , 
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so it can be expressed as tr(y4'B„_ii?^_^), where A' is 
given by the expression in the curly brackets, that pre- 
cisely gives the stated recursion formula of Eq. [241 This 
proves the proposition. □ 



Next, one turns to the detection theorem itself. This 
theorem states that any NPT entangled state will be de- 
tected in some step of the averages iteration. 

Theorem IV. 1 (Detection for the averaged iteration). 

For any detectable operator ^ with tr(p'") = 1, and 
any strictly positive operator Bi > 0, the averaged itera- 
tion process will always detect the state. 

Proof. In order to prove detection of an arbitrary state 
^ 0, it is sufhcient to show that the operators A„, 
which determine the average expectation value via Prop. 
IIV.!) become negative definite at some point in the se- 
quence. Since all operators An will necessarily commute 
with Ai = , one can easily identify the corresponding 
eigenvectors and just needs to examine the behavior of 
the corresponding eigenvalues under the iteration. As- 
sume the following spectral decompositions of the oper- 
ators 

/ = E^^I^'^X^^I+E^^I^^-X^^I' (28) 

An = J2 1^^) ("^1 + E ^1"' 1^^) ("^-1 ' (29) 
iei+ jei- 

in which the index set X_ labels all the negative eigen- 
values of , and the strictly positive eigenvalues. 

In order to show that all eigenvalues become negative 
at some point in the iteration, we proceed as follows: 
First we show that all eigenvalues on the negative sub- 
space are decreasing exponentially with the number of 
iterations n, whereas eigenvalues on the positive semidef- 
inite subspace can only increase linearly with n. This al- 
ready guarantees exponential fast divergence of tr(A„). 
Because of that, all eigenvalues on the positive subspace 
will decrease exponentially fast as well, since the factor 
tr(A„) enters in the iteration formula, cf. Eq. [24l Com- 
bined with the fact that i?i > this proves the claim. In 
detail, we show by induction the following bounds: 



«!"^ < |A+| + (n-l) 



< -|A- 




n-l 



(30) 
(31) 



for all n > 1. The induction start is trivial and we care 
about the induction step n i-^ n + I only. Let us con- 
sider the positive subspace first. In total we obtain the 



following sequence of inequalities 
' ' \ / d 



< 



< 



d 



\K\ + {n-l) 



1 - 



+ 12 



2|A+|\ , |A;| 



< |A+| + [(n + l)-l]J 



y + |2 



d 

(32) 



In the first line one has used the recursion formula of the 
operators An and the first inequality stems from the con- 
dition tr(A„) < tr(Ai) = 1, that comes from the generic 
properties of the iteration process. Prop. III. II applied to 
the average iteration process with Bi = \. In the next 
step, the induction hypothesis is employed and one ob- 
tains the final result if one upper bounds the term in 
the parenthesis. For the negative subspace, one similarly 
obtains 



< 6f ^ 1 



< 1 



2iA7n 


\ tr(A„)| 


d , 


' ^ d 


2|A7l^ 




d , 




2|A7|^ 




d 





< -IA7I 1 



|A7l 



(n+l)-l 



(33) 



Again, the recursion formula and tr(yl„) < 1 were cm- 
ployed first. The second inequality originates from the in- 
duction hypothesis, since it allows to infer |A~| < —b^f^. 
The induction step finishes with another application of 
the induction hypothesis for the last inequality. This 
concludes the proof of the theorem. □ 



V. EXAMPLES 

In order to visualize the effect of the iteration process — 
and to provide a qualitative picture of the "curvature" of 
the corresponding nonlinear entanglement witnesses — a 
simple two qubit example is sufficient. 



A. Optimized Iteration 

To investigate the optimized iteration, we consider a 
starting witness that is a convex combination between an 
optimal entanglement witness and the identity operator; 
the parameter e describes the corresponding mixedncss of 
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FIG. 2: Optimized iteration: The set of physical states is 
given by the triangular shaped region connecting the extreme 
points (0,1) and (1,0) with the origin, while all separable 
states are described by the region inside. The black triangu- 
lar, with coordinates (0.6,0.2), symbolizes the target state. 
For the last case, e — 1, each nonlinear entanglement follows 
the first linear entanglement witness (case n = 1) and "comes 
back" on the other side. 



these two operators. This class of entanglement witnesses 
with the corresponding operators Bi is given by 

W{e) = [b^bIY = e (|V.-) (^-|)^ + (1 - £)1, (34) 

Bi =i [VTT3^\i'-){^j-\ + VT^{i~yj-){^-\)] , 

with li;-) = (|01) - 1 10)) /V2. In the following we con- 
sider three different cases, the optimal entanglement wit- 
ness with e = 1, and two slightly weaker witnesses with 
e = 0.8 and e — 0.5 for different kinds of two qubit states. 

First we investigate the nonlinear improvements of 
these witnesses for a particular family of Bell-diagonal 
states given by 

p=pi (0+ 1 +P2 1 + (1 - Pi - P2) ^, (35) 

with the abbreviation |(?!)±) = (|00) ± |11)) /\/2 given m 
a standard product basis, and 1/4 denotes the totally 
mixed state. Of course, only certain parameter pairs 
actually correspond to physical states, since the corre- 
sponding operator must be positive semidefinite in order 
to form a valid density operator. In the Fig. [21 the set 
of physical states corresponds to the triangular shaped 
region. The convex set of separable states is determined 
by the partial transpose and is given by the subset within 
the set of physical states. 

The examples for the optimized iteration process are 
shown in Fig. [21 in which the target state is marked by 
the black triangular. Note that this target state is on the 
"wrong side" of the state region, i.e., it is far from being 
detected by one of the linear witness given by Eq. [34l 
As one can see, the given target state is already detected 
after a few iterations. The improvement, which char- 
acterizes the extend to which the entanglement witness 
can be "bent" over to the other side of the state region, 
depends on the original strength of the starting witness. 

Starting with a rather weak witness, the first example 
with e = 0.5, already allows to curve the witness in such a 
way that nearly all states on the other side are detected, 
while for the slightly stronger witness, the second case 
with e = 0.8, the first improvement is not yet enough to 
verify entanglement for the target state. For the optimal 
entanglement witness, the last case with £ = 1, this task 
becomes the hardest, however even in this case the tar- 
get state is detected already after the third iteration. For 
completeness, notice that starting with the identity op- 
erator, e = 0, enables detection of any given target state 
already after the first iteration. In addition note that if 
one follows the optimized iteration process, the chosen 
target states is detected, but in general it will not be 
the case that one ends up with an entanglement witness 
which is capable to witness all possible states at once. 

In a second example we investigate the number of re- 
quired iterations to detect a given target state. To this 
aim, we generated randomly (in Hilbert Schmidt norm) 
entangled two-qubit states, and computed how many it- 
erations are necessary, until they are detected starting 
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FIG. 3: Optimized iteration: The probability of detecting a 
randomly chosen entangled state in the n-th step is shown, 
for the three values of e. See text for further details. 
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FIG. 4: Random iteration: The probability of detecting a 
randomly chosen entangled state in certain intervals of iter- 
ation steps is shown, for the three values of e. See text for 
further details. 
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n=20 
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from each of the three witnesses given by Eq. [Mj The 
results are shown in Fig. [31 One can clearly see that most 
of the states are detected after less than five iterations. 



B. Random Iteration 

Let us now investigate the random iteration. For this 
aim, we consider the three witnesses given by Eq. [Ml and 
consider again randomly generated entangled states, and 
the number of iterations which are required for an itera- 
tion. The results arc plotted in Fig. [H On can sec that 
most states are detected after hundred iterations by ran- 
domly chosen unitaries. In fact, we were unable to find 
an example of a state (also for other witnesses), which is 
not detected after five hundred iterations. This suggests 
that maybe also the random iteration is complete, in the 
sense that each entangled state is detected after a finite 
number of steps. 




FIG. 5: Averaged iteration: The set of physical states and 
the set of separable states are given as in optimized iteration, 
cf. Fig. [2l In the last case, e = 1, the detected states are 
again the ones which are already detected by the first, linear 
witness (case n — 1), plus the states in the right region if the 
witness curves back. 
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C. Averaged Iteration 

Finally, Fig. [5] shows similar examples for the aver- 
aged iteration process, in which the same starting wit- 
nesses are employed as for the optimized sequence. As 
one can see, the region of detected states increases with 
the number of considered iterations, however much more 
iterations are needed in comparison to optimized itera- 
tion process. For the first example, e = 0, the identity 
operator is used as the starting operator. Using the aver- 
aged iteration method, the part of detected states is sym- 
metric with respect to the upper and lower state region, 
which shows that the corresponding witness is improved 
equally in all possible directions. This symmetry breaks 
of course if one uses one of the asymmetric entanglement 
witnesses to start with. 



VI. EXTENSIONS 

Although the iteration process and its corresponding 
results were solely discussed for the PPT criterion, the 
method directly generalizes to other trace-preserving, 
positive but not completely positive maps A; for exam- 
ple to the reduction criterion and various extensions of it 
m, , or the Choi-map [sll . All entangled states which 
violate the corresponding condition = id (8) A(p) > 
can be detected by the iteration method. In order to 
obtain the entanglement witness in the usual sense, one 
employs the adjoint map A^, defined by the property that 

tr(A(X) y) = tr(X At(y)), \JX,Y. (36) 

Therefore the entanglement witness becomes W \d® 
A'^{BB'^), which represents the general connection be- 
tween the entanglement witness and the iteration opera- 
tor B. 

Finally, due to the isomorphisms studied by de Pillis, 
Jamiolkowski, and Choi [H, H^, , we can find for a 
given witness always a positive, but not completely posi- 
tive map A, and a state such that the witness can be 
written as 

M^ = id®At(|V)(V'|), (37) 

see Ref. [UlllH for details. This is the same structure as 
the witness in Eq. [3l hence all the results can be applied. 

VII. CONCLUSION 

In conclusion, we provided a sequence of nonlinear en- 
tanglement witnesses, defined as functionals on the set of 



quantum states, which necessarily are non-negative for all 
separable states. 

Two particular iterations were investigated in more de- 
tail. In the optimized iteration one tries to optimize the 
improvement, given by the nonlinear term, according to 
a given preselected target state. Form the geometric pic- 
ture of entanglement witness, cf. Fig. [Tl this process 
corresponds to the task that one likes to "curve" a given 
witness along a certain direction. By contrast, the av- 
eraged iteration deals with the exact opposite case; one 
tries to improve the witness along all possible directions. 
The main result of the manuscript is that both iteration 
methods are successful: Any entangled state, detectable 
by the corresponding positive but not completely positive 
map, is also detected by the corresponding sequence of 
nonlinear entanglement witness, as long as the starting 
witness does not act onto a completely different subspace. 

There are several open questions which deserve further 
study. First, it remains open whether a similar detection 
statement holds for the random iteration process, as sug- 
gested by the numerical example of Sec. IV Bl Second, it 
might be desirable to extend the iteration idea beyond a 
certain positive but not completely positive map; how- 
ever this seems to require some extra knowledge about 
the structure of the quantum states applied to some par- 
ticular map. In addition, it is tempting to ask whether 
similar ideas translate to other known entanglement cri- 
teria, which do not rely on a positive but not completely 
positive maps, but which can be considered in an entan- 
glement witness form, e.o. the computable cross norm or 
realignment criterion Finally, it remains to clarify 

the possible connection between the averaged iteration 
process and the spectrum estimation idea from Ref. [13] . 
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